REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMB  No.  0704-0188 


The  public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources, 
gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of 
information,  including  suggestions  for  reducing  the  burden,  to  Department  of  Defense,  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports  (0704-0188), 
1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington,  VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any 
penalty  for  failing  to  comply  with  a  collection  of  information  if  it  does  not  display  a  currently  valid  OMB  control  number. 

PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 


1.  REPORT  DATE  (DD-MM-YYYY)  2.  REPORT  TYPE 

06/29/2006  Final 


4.  TITLE  AND  SUBTITLE 

Coherent  Structures  and  Chaos  Control  in  High-Power  Microwave  Devices 


3.  DATES  COVERED  (From  -  To) 

1  April  2003-  30  March  2006 


5a.  CONTRACT  NUMBER 


5b.  GRANT  NUMBER 

F49620-03- 1-0230 


5c.  PROGRAM  ELEMENT  NUMBER 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Massachusetts  Institute  of  Technology 
Cambridge,  MA  02139 
Kensington,  MD  20895 


5d.  PROJECT  NUMBER 


5e.  TASK  NUMBER 


5f.  WORK  UNIT  NUMBER 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

Air  Force  Office  of  Scientific  Research 
4015  Wilson  Blvd 
Mail  Room  713 
Arlington,  VA  22203 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Distribution  A;  distribution  unlimited 


10.  SPONSOR/MONITOR  S  ACRONYM(S) 

AFOSR 


AFRL-SR-AR-TR-06-0290 


13.  SUPPLEMENTARY  NOTES 


14.  ABSTRACT 

This  report  summarizes  research  results  obtained  under  the  auspices  of  Air  Force  of  Scientific  Research,  Grant  No. 

F49620-03- 1-0230  (Chen,  2003).  In  particular,  we  conducted  vigorous  theoretical  and  computational  investigations  of  coherent 
structures  and  chaos  in  a  wide  range  of  intense  electron  beam  plasmas  relevant  to  the  development 
of  high-power  microwave  sources  for  directed  energy  applications.  The  following  is  a  brief  summary  of  our  research 
accomplishments  in  selected  areas,  while  detailed  findings  are  described  in  the  preprints  cited  in  this  report. 


15.  SUBJECT  TERMS 


20060727323 


16.  SECURITY  CLASSIFICATION  OF: 

a.  REPORT 

b.  ABSTRACT 

c.  THIS  PAGE 

U 

U 

u 

17.  LIMITATION  OF  18.  NUMBER  19a.  NAME  OF  RESPONSIBLE  PERSON 
ABSTRACT  OF 

PAGES  - 

TJTJ  19b.  TELEPHONE  NUMBER  (Include  area  code) 


Standard  Form  298  (Rev.  8/98) 

Prescribed  by  ANSI  Std.  Z39.18 


Final  Report 

Coherent  Structures  and  Chaos  Control  in  High-Power  Microwave  Devices 

AFOSR  Grant  No.  F49620-03-1-0230 

Submitted  to: 

Dr.  Arje  Nachman 
Program  Manager 

Air  Force  Office  of  Scientific  Research 
875  North  Randolph  Road 
Suite  325,  Room  3112 
Arlington,  VA  22203 
Phone:  703-696-8427 
FAX:  703-696-8450 

Submitted  by 

Plasma  Science  and  Fusion  Center 
Massachusetts  Institute  of  Technology 
Cambridge,  Massachusetts  02139 

Principal  Investigator 
Dr.  Chiping  Chen 


Grant  Period 

April  1,  2003  -  March  30,  2006 


Table  of  Contents 


Cover  Page  1 

Table  of  Contents  2 

Final  Report  .  3 

1 .  High-Power  Magnetron  Research  3 


1.1.  Discovery  of  Vortex  Structures  in  Magnetrons 

1.1.1.  Guiding-Center  Theory 

1.1.2.  Cold-Fluid  Theory 

1.2.  Formulation  of  Small-Signal  Gain  Theory  of  Multiresonator 
Cylindrical  Magnetrons 

2.  High-Power  Klystron  Research  9 

2.1.  Determination  of  the  Current  Limit  on  the  Confinement  of  Finite-Size 
Bunched  Pencil  Beams  in  High-Power  Relativistic  Klystrons 

2.2.  Exploration  of  the  Possibility  of  Magnetic  Cusp  Formation  in  Highly 
Bunched  Annular  Beams  in  High-Power  Relativistic  Klystrons 

3.  Development  of  Ellipse-Shaped  Ribbon-Beam  Theory  for  HPM  Device 

Applications  12 

3.1.  Theory  of  Elliptic  Beam  Formation 

3.2.  Theory  of  Elliptic  Beam  Focusing  and  Transport 

3.3.  Spin-off  of  Ribbon-Beam  Technology  for  Commercial,  Defense  and 
Scientific  Applications 

4.  Chaotic  Particle  Motion  and  Halo  Formation  in  Small-Aperture  Systems  20 

5.  Research  Initiative  on  Relativistic  Electron  Flows  and  Diodes  for  HPM  Device 

Applications  '22 

6.  References  23 


2 


Final  Report 

Coherent  Structures  and  Chaos  Control  in  High-Power  Microwave  Devices 
AFOSR  Grant  No.  F49620-03-1-0230 

This  report  summarizes  research  results  obtained  under  the  auspices  of  Air  Force  of 
Scientific  Research,  Grant  No.  F49620-03-1-0230  (Chen,  2003).  In  particular,  we 
conducted  vigorous  theoretical  and  computational  investigations  of  coherent  structures 
and  chaos  in  a  wide  range  of  intense  electron  beam  plasmas  relevant  to  the  development 
of  high-power  microwave  sources  for  directed  energy  applications. 

The  following  is  a  brief  summary  of  our  research  accomplishments  in  selected  areas, 
while  detailed  findings  are  described  in  the  preprints  cited  in  this  report. 

1.  High-Power  Magnetron  Research 

Under  the  auspices  of  the  present  grant  (Chen,  2003),  we  developed  equilibrium  theories 
of  electron  flow  in  magnetrons  and  crossed  field  devices,  taking  into  account  the 
influence  of  the  corrugation  on  the  anode  on  the  electron  flow  properties.  Under  certain 
conditions,  both  theories  predicted  the  existence  of  vortex  structures  induced  by  the 
corrugation  on  the  anode  (Davies  and  Chen,  2006;  Bhatt  and  Chen,  2003  and  2004). 
Vortex  formation  is  a  new  phenomenon  which  had  not  been  discussed  in  the  literature  in 
the  context  of  magnetrons  and  crossed  field  devices.  Since  the  report  of  vortex  formation 
(Chen,  2004)  at  Magnetron/PIC  Simulation  Workshop  held  at  AFRL,  May  2-3,  2004, 
several  researchers  at  AFRL  (Cartwright,  2005)  and  elsewhere  (Bosman,  et  al.,  2005)  are 
exploring  use  of  non-axisymmetric  cathodes  to  improve  magnetron  performance. 

We  also  formulated  a  small-signal  gain  theory  for  magnetrons  (Chen,  2005),  using 
both  the  Floquet  theorem  and  the  Maxwell  equations.  Once  it  is  fully  developed,  it  will 
advance  the  existing  magnetron  gain  theories  (Riyopoulos,  1998;  Kaup,  2005,  and 
references  therein)  which  do  not  include  them  simultaneously. 

1.1.  Discovery  of  Vortex  Structures  in  Magnetrons  (Davies  and  Chen,  2006;  Bhatt 
and  Chen,  2003  and  2004) 

Until  our  work  (Davies  and  Chen,  2006,  Bhatt  and  Chen,  2003  and  2005),  the  influence 
of  the  anode  corrugation  on  the  equilibrium  electron  flow  properties  had  not  been 
discussed  in  the  magnetron  literature.  In  fact,  all  conventional  treatments  of  the 
equilibrium  electron  flows  in  magnetrons  had  ignored  effects  of  the  corrugations  at  the 
anode  (see,  for  example,  Davidson,  et  al.,  1991). 

We  developed  two  equilibrium  theories  to  describe  the  equilibrium  electron  flow  in 
magnetrons  and  crossed  field  devices.  One  of  the  equilibrium  theories  was  based  on  a 
guiding-center  model  in  which  the  electron  flow  velocity  is  given  by  the  ExB  drift 
velocity  (Davies  and  Chen,  2006).  The  other  was  a  general  treatment  in  which  the 
electron  flow  velocity  is  determined  from  the  cold-fluid  force  equation  (Bhatt  and  Chen, 
2003  and  2004).  Both  theories  predicted  the  formation  of  periodic  vortex  structures  in  the 
equilibrium  electron  flow. 

By  definition,  an  equilibrium  flow  corresponds  to  a  state  in  magnetrons  and  crossed- 
field  devices  in  the  absence  of  rf  oscillations.  While  it  is  an  ideal  situation  in  an  actual 
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device,  it  is  of  critical  importance  in  order  for  us  to  develop  a  better  understanding  of 
magnetrons  and  crossed-field  devices. 

1.1.1.  Guiding-Center  Theory  (Davies  and  Chen,  2006) 

We  have  made  an  analysis  of  non-axisymmetric  ExB  equilibrium  flow  in  a  crossed- 
field  apparatus  (Davies  and  Chen,  2006).  The  apparatus  consists  of  a  conducting  circular- 
cylindrical  cathode  encompassed  by  either  a  corrugated  conducting  anode  or  by  a 
periodically  biased,  circular-cylindrical  wall.  A  uniform  magnetic  field  is  directed 
parallel  to  the  cathode  axis.  With  the  aid  of  perturbation  theory,  the  guiding-center 
approximation,  and  the  assumption  of  a  constant  electron  density  in  the  electron  layer,  we 
determine  equipotential  surfaces,  electron  trajectories,  and  the  corrugated  boundary  of  the 
electron  layer. 

An  interesting  result  is  the  appearance  of  vortex  structures  in  regions  of  negative 
potential  near  the  cathode  surface.  The  scaling  properties  of  the  width  and  potential 
depth  of  the  vortex  are  studied  for  the  small-amplitude  corrugation  regime.  For 
sufficiently  thin  electron  layers,  results  are  shown  to  be  applicable  to  relatively  high 
density  regimes  as  well  as  to  low  density  regimes. 

While  detailed  results  of  this  research  are  described  in  the  manuscript  “Vortex 
structure  in  a  high-density  non-axisymmetric  ExB  equilibrium  flow”  by  J.A.  Davies 
and  C.  Chen  at  the  appendix  of  this  report,  we  highlight  below  some  the  key  features  of 
the  vortex  structures. 

Figure  1.1  shows  a  schematic  of  a  magnetron  in  a  two-dimensional  configuration, 
where  the  corrugated  anode  corresponds  to  an  equipotential  surface.  Rather  than 
designing  the  anode  corrugations  to  achieve  desired  electromagnetic  properties  in  an 
actual  device  in  the  conventional  design  practice,  we  determine  the  corrugated  anode 
self-consistently  from  the  theory,  as  an  inverse  problem.  Because  there  is  considerable 
freedom  for  us  to  choose  the  anode  within  our  theory,  we  can  fine-tune  the  corrugated 
anode  give  us  the  desired  electromagnetic  properties  in  magnetron  design. 

Figure  1.2  shows  the  equilibrium  electron  flow  for  the  case  of  a  uniform  density 
profile  extending  from  the  cathode  to  the  periodic  electron  layer  boundary.  In  Fig.  1.2, 
the  normalization  is  chosen  as  p-r!a  and  C,  =0/4moea2 ,  where  a  is  the  cathode 
radius,  -e  is  the  electron  charge,  <f)  is  the  electrostatic  potential,  and  n0  is  the  electron 
density.  The  cathode  is  at  p- 1  =  0,  whereas  the  anode  is  at  £  =  =  0.03 .  The  solid 

curve  with  ^  =  0.00301  represents  the  electron  layer  boundary.  Electrons  move  along  the 
equipotential  curves.  The  close-up  shows  one  of  6  vortices  near  the  cathode. 

In  general,  vortex  structures  form  near  the  space-charge-limited  emission  condition, 
which,  in  our  theory,  corresponds  to  the  condition  that  the  radial  electric  field  averaged 
over  the  azimuthal  angle  vanishes  at  the  cathode. 

We  have  also  shown  that  vortex  structures  occur  when  the  electron  density  is  non- 
uniform,  which  will  be  a  subject  of  our  second  paper  on  vortex  structures  in  the 
equilibrium  electron  flow  in  magnetrons  and  crossed-field  devices. 
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Figure  1.1:  Schematic  of  a  magnetron  in  a  two-dimensional  configuration.  The 
cathode  is  located  at  r  =  a,  whereas  the  corrugated  anode  is  at  r  =  rAN{0).  The 
electron  layer,  which  extends  from  the  cathode  to  the  e-layer  boundary,  is 
confined  by  the  applied  magnetic  field  50ez  (from  Davies  and  Chen,  2006). 


Figure  1.2:  Equipotential  curves  in  the  equilibrium  electron  flow  (top)  and  their 
close-up  (bottom).  The  cathode  is  at  p- 1  =  0,  whereas  the  anode  is  at 
^  =  0.03  .  The  solid  curve  with  C,  =  0.00301  represents  the  electron  layer 

boundary.  Electrons  move  along  the  equipotential  curves.  The  close-up  shows  one 
of  6  vortices  near  the  cathode  (from  Davies  and  Chen,  2006). 
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1.1.2.  Cold-Fluid  Theory  (Bhatt  and  Chen,  2003  and  2004) 


The  guiding-center  model  discussed  in  Sec.  1.1.1  is  applicable  to  non-relativistic 
magnetrons  and  crossed-field  devices  with  thin  electron  layers,  i.e.,  r/a- 1«1.  We 
also  investigated  vortex  structures  in  the  equilibrium  using  a  cold-fluid  model.  In  the 
cold-fluid  model,  the  restriction  rt a- 1 «  1  was  removed. 

Because  the  electron  flow  velocity  is  given  by  the  cold-fluid  force  balance  equation 
rather  than  the  simple  E  x  B  drift  (which  is  incompressible  flow,  i.e.,  divergence-free 
flow),  the  electron  motion  is  considerably  more  complicated,  in  general.  The  problem  is 
still  tractable. 

In  particular,  we  solved  the  2D  cold-fluid  equations  by  introducing  a  generalized 
boundary-conformal  coordinate  system  and  scalars  representing  curl-free  and  divergence- 
free  velocity  fields.  The  cold-fluid  equations  are: 

0  =  {o)c+  V2W)V20  +  V(V2W)-  VO  +  e2  •  [v'F  x  v(v2vf)]  (1.1) 

and 

0  =  (v20  +  V  •  V)[I  V2 V2  -  [coc  +  V2w)v2W  -  v(v2vf)-  VW  -  e2  •  [v(v2vp)x  Vo]j  (1 .2) 

where  V  =  VO  +  ez  x  VW  and  coc  =  qB/mc . 

For  a  uniform  magnetic  field  and  uniform  energy  profile  across  the  fluid,  we  found 
that  V2W  =  -coc  =  const  ,  which  corresponds  to  a  generalization  of  Brilliouin  flow 
(Davidson,  1990).  In  this  situation,  the  2D  cold-fluid  equations  reduce  to  a  single 
nonlinear  partial  differential  equation,  i.e., 

0  -  [v20  +  (VO  +  e2  x  VT')  •  v]v2  [(w)2  +  (VO)2  +  2e2  •  (VW  x  VO)J.  (1.3) 

For  small  departures  of  the  electron  flow  from  azimuthal  symmetry,  we  solved  the 
nonlinear  equation  (1.3)  perturbatively. 

As  an  example,  we  considered  a  simple  case  where  the  inner  boundary  of  the  electron 
layer  is  azimuthally  symmetric.  In  terms  of  units  normalized  to  this  inner  radius,  we 
represented  this  surface  by  r  =  a{0)  =  1 .  On  the  outer  boundary  of  the  electron  layer,  we 
introduced  a  small  ripple,  describing  it  by  r  -  b($)=  1.8  x  [l  +  0.01  cos(2f?)].  Solving  for 
the  fluid  flow  in  this  thick  annulus,  we  found  areas  where  the  velocity  field  reverses  on 
itself,  forming  fluid  vortices  as  shown  in  Fig.  1.3.  These  vortices  are  also  regions  of 
potential  minima,  providing  opportunities  for  ion  trapping  and  the  associated  noise 
mechanisms  (Yamamoto,  et  al.,  1987). 

Our  discovery  of  vortex  structures  in  magnetrons  and  crossed-field  amplifiers 
provides  new  insight  into  a  number  of  areas  in  magnetron  and  CFA  research,  including 

a)  Elimination  of  ion  trapping  and  ion  noise, 

b)  Elimination  of  turbulence  and  noise, 

c)  Improvement  of  magnetron  efficiency  (Cartwright,  2005),  and 
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Figure  1.3:  Vector  field  plots  of  the  flow  velocity  as  computed  from  the 
perturbation  theory  for  a  small  ripple  on  the  outer  fluid  boundary  surface.  Color 
indicates  the  magnitude  of  the  fluid  velocity  at  a  given  location.  Shown  on  the  top 
is  the  entire  annulus.  A  close  up  of  the  rectangular  region  on  the  annulus  (bottom) 
reveals  a  portion  of  a  vortex  centered  at  jc  =  1.06  and  y  =  0  (from  Bhatt  and 
Chen,  2003). 
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d)  Suppression  of  mode  competition  (Bosman,  et  al,  2005). 
Some  of  these  items  will  be  the  subject  of  our  future  investigations. 


1.2.  Formulation  of  Small-Signal  Gain  Theory  of  Multiresonator  Cylindrical 
Magnetrons  (Chen,  2005) 

We  formulated  a  small-signal  gain  theory  of  magnetrons  for  a  muliresonator  cylindrical 
magnetron  (Fig.  1.4)  using  the  linearized  cold-fluid  equations,  the  full  Maxwell  equations, 
and  the  Floquet  theorem.  The  gain  equation  is  expressed  as  an  eigenvalue  problem 
defined  by  (Chen,  2005) 


1  d 
r  dr 


r 

1  -  co2r2  /  c2{l  +  kN)2 


(l  +  Xrk,) 


dM>kl 

dr 


(l  +  kN)2 


0  +  Xeki)#&ki  =  0 


(1.4) 


under  appropriate  boundary  conditions.  In  Eq.  (1.4),  dd>k,  is  proportional  to  the  electric 
field  perturbation,  Xm  and  are  susceptibility  functions,  a>  is  the  frequency  of  the 
perturbation,  N  is  the  number  of  vanes  on  the  anode,  the  index  /  with  /  =  1,  2 N 
designates  the  azimuthal  mode  number,  and  the  index  k  with  £  =  0,±1,...  designates 
the  Floquet  mode  number. 

In  principle,  solving  the  eigenvalue  problem  will  allow  us  to  determine  the  small- 
signal  gain  in  a  magnetron  for  any  given  equilibrium  electron  flow  profile.  As  we  will 
demonstrate  in  Sec.  3.1,  solving  the  eigenvalue  equation  in  the  absence  of  the  electron 
layer  yields  a  vacuum  dispersion  relation  which  agrees  with  experimental  measurements. 
We  plan  to  solve  the  eigenvalue  problem  with  the  presence  of  the  electron  layer  and 
extend  the  small-gain  theory  to  the  relativistic  regime. 


Figure  1.4:  Schematic  of  a  multiresonator  cylindrical  magnetron  in  a  two- 
dimensional  configuration.  The  cathode  is  located  at  r  =  a,  whereas  the  anode  is 
a  corrugated  slow-wave  structure.  The  electron  layer  is  confined  by  the  applied 
magnetic  field  B0ez  (from  Chen,  2005). 
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2.  High-Power  Klystron  Research 

In  the  operation  of  high-power  relativistic  klystrons,  it  is  critical  to  confine  the  electron 
beam  radially  as  it  becomes  sharply  bunched  along  the  direction  of  propagation. 
Typically,  a  high-power  relativistic  klystron  employs  either  a  pencil  relativistic  electron 
beam  or  an  annular  relativistic  electron  beam. 

Under  the  auspices  of  the  grant  (Chen,  2003a),  we  extended  our  treatment  of  the 
bunched  beam  confinement  from  our  previous  point-charge  model  (Hess  and  Chen,  2000 
and  2002b,  Hess,  2002)  to  a  finite-size  bunched  beam  model  (Hess  and  Chen,  2004),  and 
successfully  applied  it  to  improving  our  understanding  of  the  operational  limit  of  high- 
power  relativistic  klystrons  (Hess  and  Chen,  2004).  We  also  explored  the  possibility  of 
magnetic  cusp  formation  in  highly  bunched  annular  beams  (Hess  and  Chen,  2006). 

2.1.  Determination  of  the  Current  Limit  on  the  Confinement  of  Finite-Size  Bunched 
Pencil  Beams  in  High-Power  Relativistic  Klystrons  (Hess  and  Chen,  2004) 

We  derived  the  space-charge  limit  for  a  train  of  periodic  finite-size  charged  bunches 
which  are  slightly  displaced  from  the  axis  of  a  circular,  perfectly  conducting  pipe.  The 
general  calculation  assumed  that  the  bunches,  while  periodic,  may  have  arbitrary  radial 
and  longitudinal  density  distributions  with  respect  to  the  beam  axis.  We  obtained  results 
for  a  few  specific  distributions,  as  well  as  the  limits  of  unbunched  and  single  bunch 
beams.  Finally,  we  showed  how  the  bunched  beam  space-charge  limit  along  with  the 
radial  rms  envelope  equation  for  a  bunched  beam  is  applied  to  predict  operational  limits 
of  periodic  permanent  magnet  (PPM)  focusing  klystrons  which  were  designed  and  tested 
for  Linear  Collider  (LC). 

The  results  are  summarized  in  Table  2.1  (from  Hess  and  Chen,  2004).  Table  2.1  lists 
the  key  parameters  for  the  five  PPM  klystrons  (Sprehn,  et  al,  1998;  Sprehn,  et  al,  2000; 
Matsumoto,  et  al,  2001;  Chin,  et  al,  2001;  SLAC  Website  http://www- 
group.slac.stanford.edu/kly/ppm_klystron/ppm_kly.html):  operating  frequency  /  , 
average  beam  current  lh ,  beam  energy  yh ,  the  rms  magnetic  field  strength  Brms ,  and  the 
pipe  radius  a.  Table  2.1  also  shows  the  values  for  a  =  2mjybL  and  the  experimental 
space  charge  (current)  parameter  s  ,  which  can  be  computed  using  the  formulas 
L  =  Pbc / /  and  Nh  =  1,,/ef .  The  theoretical  limit  is  se  /se(0,0).  If  the  experiment  value 
seexp/se(0,0)  is  greater  than  se  /^(O.O)  in  a  klystron,  then  substantial  beam  loss  occurs. 
If  the  experimental  value  s  /  se  (0,0)  is  less  than  5c/^(0,0)  in  a  klystron,  then  small 

beam  loss  could  still  occurs  due  to  beam  halo  formation  (Chen  and  Pakter,  2000;  Pakter 
and  Chen,  2000).  Our  results  show  good  agreement  with  the  experimental  observations. 

2.2.  Exploration  of  the  Possibility  of  Magnetic  Cusp  in  Highly  Bunched  Annular 
Electron  Beams  in  High-Power  Klystrons  (Hess  and  Chen,  2006) 

We  explored  the  possibility  of  magnetic  cusp  formation  in  high-power  klystrons  driven 
by  high-intensity  annular  beams,  such  as  the  relativistic  klystron  oscillator  which  is 
currently  being  studied  experimentally  at  Air  Force  Research  Laboratory. 
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For  a  thin  annular  disk,  the  equilibrium  force  balance  equation  has  the  dimensionless 
form  (Hess  and  Chen,  2006) 


/u^\~r2(o7 


/t  °° 

=  5  -  -  £  J,  (y0m^)coth[^0m  /a] 

-  Aa>\ £ Km J 0 (jQJ )coth [nj0m /a]  _  — 

lm=l  A  m=l  Jom 


2  a2  fyyyil/’oJ^oK)] 


(2.1) 


where  ct){r)=V™'(r)/r  ,  a>  =0x1/0  ,  a  =  ana2  ,  K  =  Kna2/c  ,  n  =  cdcajc  , 
A  =  2Nbe2  / meca>ca2 ,  and  coc  =eBalmec .  Equation  (2.2.1)  predicted  the  formation  of  a 
magnet  cusp  as  the  beam  current  is  increased,  as  shown  in  Fig.  2.2.1. 

Figure  2.1  shows  plots  of  acolc  and  Bself  / B0  vs.  r/o  for  //  =  10.7  and  two  cases  of 
A  =  0.05  and  A  =  0.1 1 15 .  The  surface  number  density  profile  is  chosen  as 

0,  r<r„ 

*r(r)  =  ■  3  (ra  -  r)(r  -  rt )/  nFS2 ,  rt<r<r0,  (2.2) 

l  0=  ra<r, 

where  r  =  {ra  +  rt)/2  =  0.8a  is  an  average  beam  radius  and  S  =  r0—ri  =  0.2a  is  the  beam 
width.  Notice  that  as  A  increases  to  a  critical  value  A  =  0.1 1 16,  Bse,f  (r)/B0  forms  a 
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Figure  2.1:  Plots  of  (top)  aco  /  c  and  Bsef  /B0  vs.  r  I  a  predicted  by  Eq.  (2.1)  for  //  =  10.7 
and  two  cases  of  X  =  0.05  and  X  =  0.1115,  using  the  surface  density  profile  defined  in 
Eq.  (2.2)  (from  Hess  and  Chen,  2006). 


cusp  shape  at  a  specific  point  r*  =  0.7881,  with  Bsef  (r  *)/B0  =  -1 .  Near  the  cusp,  we 
found  numerically  that  the  magnetic  field  scales  approximately  as 
B'elf  (r)+  B0  oc  (r*  -rf'  with  ax  =  0.5  for  r  </  ,  and  Bself(r)+B0  °c  (r-r*)  2  with 
a2  s  0.73  for  r>r\  Beyond  the  cusp  with  X  >  Xcmp ,  we  were  unable  to  find  numerical 

solutions  to  the  equilibrium  momentum  equation. 

The  physics  of  magnetic  cusp  formation  is  that  as  the  beam  current  increases,  the 
rotation  of  the  electron  flow  becomes  relativistic  and  produces  a  diamagnetic  field  which 
conceals  the  applied  magnetic  field. 

While  magnetic  cusp  formation  is  a  new  phenomenon  predicted  by  our  theory,  it 
remains  to  be  seen  whether  such  phenomena  can  be  observed  in  annular-beam  driven 
klystron  experiments.  A  key  question  would  be  whether  the  electrons  will  be  sufficiently 
bunched  along  the  direction  of  beam  propagation,  which  will  be  a  subject  of  our  future 
investigations. 
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3.  Development  of  Ellipse-Shaped  Ribbon-Beam  Theory  for  High-Power  Microwave 
Device  Applications 

Existing  HPM  devices  employ  either  a  pencil  electron  beam  or  an  annular  electron  beam. 
One  of  the  main  objectives  in  high-power  microwave  (HPM)  research  is  to  develop 
innovative  science  and  technology  which  could  lead  to  improve  the  efficiency,  output 
power  level,  and  other  performance  (such  as  weight  reduction)  of  high-power  microwave 
(HPM)  devices. 

A  promising  approach  is  to  use  a  large-aspect-ratio  ellipse-shaped  electron  beam 
rather  than  the  conventional  pencil  or  annular  electron  beam  as  the  energy  source. 
Compared  with  a  conventional  HPM  device,  a  HPM  device  powered  by  an  ellipse-shaped 
electron  beam  has  the  following  attractive  features: 

a)  Power  scales  as  instead  of  f~2  (where  /  is  the  frequency),  thus  more 
power; 

b)  Low  effective  beam  perveance,  thus  higher  efficiency; 

c)  Use  of  permanent  magnets  for  beam  focusing,  thus  lower  energy  consumption; 
and 

d)  Low  cathode  loading,  thus  longer  device  lifetime. 

However,  there  was  lack  of  a  fundamental  understanding  of  the  generation  (or  formation) 
of  high-power  ellipse-shape  electron  beam.  Indeed,  few  analytic  results  were  available 
until  the  breakthrough  we  made  recently  in  the  development  of  three-dimensional  theory 
and  three-dimensional  simulation  of  the  generation  (or  formation)  of  ellipse-shaped 
ribbon  electron  beams. 

The  formation  and  transport  of  ribbon  (sheet)  electron  beams,  which  is  three- 
dimensional  in  nature,  had  presented  major  challenges  in  the  realization  of  ribbon-beam 
microwave  devices.  Under  the  auspices  of  the  present  grant  (Chen,  2003a),  we  made 
major  breakthroughs  in  the  theoretical  understanding  of  ribbon-beam  formation  and 
transport. 

3.1.  Theory  of  Elliptic  Beam  Formation  (Bhatt  and  Chen,  2005;  Bhatt,  Bemis  and 
Chen,  2005  and  2006) 

We  published  an  important  paper  on  a  three-dimensional  theory  of  non-relativistic, 
laminar,  elliptic  beam  formation  (Bhatt  and  Chen,  2005).  The  beam  was  shown  to  have 
the  Child-Langmuir  form  (Bhatt  and  Chen,  2005).  An  analytic  expression  for  the 
electrostatic  potential  outside  the  beam  was  derived.  Equipotentials  corresponding  to 
electrode  surfaces  were  computed  numerically.  The  effectiveness  of  the  electric  field 
formed  by  the  electrodes  in  focusing  and  preserving  the  elliptic  Child-Langmuir  beam 
was  verified  via  3D  ray-tracing  simulations. 

Figure  3.1  shows  a  schematic  of  an  elliptic  beam  and  the  elliptical  coordinate  system 
used  in  our  theory.  Figure  3.2  shows  the  contours  of  the  potentials  for  the  cathode  and 
anode  focusing  electrodes.  Figure  3.3  shows  the  OMNITRAK  simulation  of  the  laminar 
flow  electron  trajectories  for  the  case  shown  in  Fig.  3.2.  The  elliptic  beam  shown  in  Figs. 
3.2  and  3.3  has  an  aspect  ratio  of  10:1. 
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Figure  3.1:  Schematic  of  an  elliptic  electron  diode  and  elliptical  coordinate  system 
(from  Bhatt  and  Chen,  2005). 


Figure  3.2:  Potential  contours  at  (a)  cathode  and  (b)  anode  focusing  electrodes 
(from  Bhatt  and  Chen,  2005). 
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Figure  3.3:  Laminar  flow  electron  trajectories  from  a  3D  OmniTrak  simulation 
projected  to  (a)  vertical  plane  (y,z )  and  (b)  horizontal  plane  {x,z)  (from  Bhatt 
and  Chen,  2005). 


The  treatment  was  extended  (Bhatt,  Bemis  and  Chen,  2005  and  2006)  to  consider 
whether  the  diode  geometry  solutions  thus  obtained  are  robust  to  perturbations  and 
limitations  of  the  sort  likely  to  be  encountered  in  a  realistic  device:  finite  extent,  part 
misalignment,  tolerances  for  mechanical  and  thermal  stresses,  etc.  Analytic  and  semi- 
analytic  estimates  were  obtained  along  with  simulations  utilizing  the  3D  trajectory  code, 
OMNITRAK.  It  was  found  that  the  elliptic-beam  solution  is  quite  stable  and  robust,  and 
its  desirable  properties  can  be  maintained  in  a  realistic  diode. 

For  his  joint  paper  “Three-Dimensional  Theory  and  Simulation  of  Large-Aspect-Ratio 
Ellipse-Shaped  Charged-Particle  Beam  Gun”  with  Tom  Bemis  and  Chiping  Chen  (2005), 
our  graduate  student  Mr.  Ronak  Bhatt  received  First-Place  Award  for  an  Outstanding 
Technical  Paper  by  Student  at  Particle  Accelerator  Conference,  May  16-20,  2005,  which 
over  1000  scientists  from  all  over  the  world  attended.  Mr.  Bhatt  also  presented  an  invited 
paper  at  ICOPS  2005,  Monterey,  CA  (Bhatt,  Bemis  and  Chen,  2005). 
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3.2.  Theory  of  Elliptic  Beam  Focusing  and  Transport  (Zhou,  Bhatt  and  Chen,  2006) 

For  HPM  devices,  the  transport  and  focusing  of  an  ellipse-shaped  electron  beam  is 
another  important  area  which  requires  research  and  development.  Under  the  auspices  of 
the  current  research,  we  showed  that  there  exists  an  exact  paraxial  cold-fluid  equilibrium 
of  a  high-intensity,  space-charge-dominated  charged-particle  beam  with  a  periodically 
twisted  elliptic  cross  section  in  a  non-axisymmetric  periodic  magnetic  field.  Generalized 
envelope  equations,  which  determine  the  beam  envelopes,  ellipse  orientation,  density, 
and  internal  flow  velocity  profiles,  were  derived.  Nonrelativistic  and  relativistic  examples 
of  such  beam  equilibria  were  discussed.  The  equilibrium  and  stability  of  such  beams  were 
demonstrated  by  self-consistent  particle-in-cell  (PIC)  simulations. 

Our  theory  (Zhou,  Bhatt  and  Chen,  2006)  describes  a  high-intensity,  space-charge- 
dominated  unbunched  beam,  in  which  kinetic  (emittance)  effects  are  negligibly  small. 
The  focusing  field  is  a  three-dimensional  (3D)  non-axisymmetric  periodic  magnetic  field 
which  can  be  described  in  the  paraxial  approximation  as 


B“'(x)  »  Bq 


v0x 


cos(k0s)xex  +  —  cos^s^ye  -sin(A:05)ez 


(3.1) 


where  k0  =  2 n/S,  klx  +  kly  =  k2 ,  and  S  is  the  axial  periodicity  length.  The  3D  magnetic 

field  in  Eq.  (2.3.1)  is  fully  specified  by  the  three  parameters  B0,  S  and  k0x/k0y.  The 

cold-fluid  equilibrium  equations  in  the  paraxial  approximation  are  (Zhou,  Bhatt  and  Chen 
2005): 


Pbcj-  +  Vx.(«,V1)  =  0, 
os 


(3.2) 


f  pbc—  +  VX  — 

^  os  j 


Vx  = 
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,(3.3) 

.(3-4) 


which  are  nonlinear. 

The  beam  can  be  adequately  described  by  cold-fluid  equations  since  the  kinetic 
effects  are  negligibly  small  in  the  beam  that  we  are  considering.  We  find  that  a  well 
confined  unbunched  beam  exists  in  such  a  focusing  system.  Inside  the  ellipse-shaped 
beam  boundary,  the  beam  has  the  uniform  density  profile 


«a(xi>s)  =  - 


0 


y 


a2(s) 


b2(s\ 


(3.5) 


and  the  linear  flow  velocity  profile 
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Figure  3.4:  Untwisted  and  twisted  coordinate  systems  (from  Zhou,  Bhatt  and 
Chen,  2006). 


Vj.(xx,  $) = \mAs)*  -  <*As)y]Pb<*s  +  \py(s)y + ay{s)x\pbch  >  (3-6) 

where  xx  =  xe~  +  ye~  is  a  transverse  displacement  in  the  twisted  coordinate  system 
illustrated  in  Fig.  3.4;  0(s)  is  the  twist  angle  of  the  ellipse;  ©(x)  =  l  if  x>0  and 
@(x)  =  0  if  x  <  0 . 

It  can  be  shown  that  both  the  equilibrium  continuity  equation  and  force  equation  are 
satisfied  if  the  dynamical  variables  0(5),  b(s) ,  /2x{s)=a~xdal  ds  ,  /jy(s)  =  b~'db/  ds  , 
0^(5),  a  (s)  and  0(s)  obey  the  generalized  beam  envelope  equations  (Zhou,  Bhatt  and 


Chen,  2006) 
dda_ 
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where  xJk~  =  qBj2ybpbmc 2  and  K  s  2q2N b/ ylPlmc2  . 
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Ta 

)le  3.1  System  parameters  for  elliptic-beam  examples 

Parameter 

Nonrelativistic 

(RBA) 

Mildly 

Relativistic 

(RBA) 

Fully 

Relativistic 

(RBK) 

Fully 

Relativistic 

(RBK) 

Application 

Wireless 

Communication 

Missile 

Defense 

International 

Linear 

Collider 

Directed 

Energy 

1.95 

2.8 

1.3 

1.3 

RF  Power 
(kW) 

500 

(pulsed) 

■ 

1 

200,000 

(pulsed) 

Current  (A) 

0.11 

18.5 

85.5 

450  ’ 

Voltage  (kV) 

2.29 

45.0 

198.5 

600 

S  (cm) 

1.912 

2.626 

2.2 

3.5 

kOx/kOy 

1.60 

1.44 

1.52 

1.51 

B0  (kG) 

0.337 

1.099 

2.40 

2.50 

a/b 

6.0 

4.0 

5.0 

5.0 

a  (cm) 

0.373 

0.585 

0.425 

0.663 

Qmax  (deg) 

10.4 

11.5 

8.8 

7.2 

A  numerical  module  in  the  MIT  2D  Periodically  Focused  Beam  (PFB2D)  code  was 
developed  to  solve  the  generalized  envelope  equations  (3.7)-(3.1 1).  The  matched  solution 
was  determined  with  Newton’s  method.  The  beam  equilibria  predicted  by  the  generalized 
envelope  equations  were  verified  by  2D  particle-in-cell  (PIC)  simulations  using  the 
PFB2D  code.  In  the  PFB2D  simulations,  we  used  the  paraxial  field  in  Eq.  (3.1),  typically 
5  x  105  particles,  a  square  grid  with  400  x  400  cells,  and  a  square  conducting  pipe  with  a 
full  width  which  is  3  times  the  semi  major  axis  of  the  beam. 

As  an  example,  we  have  considered  a  relativistic  ribbon  beam  with  voltage 
Vh  =198.5  keV  ,  current  Ih  =85.5  A,  aspect  ratio  a/b  =  5  ,  and  non-axi symmetric 

periodic  permanent  magnet  focusing  with  B0  -  2.4  kG ,  S  =  2.2  cm,  and  k0y/k0x  =  1.52 

(see  Table  3.1,  Column  4).  We  propose  to  use  it  in  a  10  MW  L-Band  ribbon-beam 
klystron  (RBK)  for  a  linear  collider  (LC). 

For  such  a  system  the  matched  solution  of  the  generalized  envelope  equations  is 
calculated  numerically  as  shown  in  Figs.  3.5(a)-(c)  with  the  corresponding  parameters: 

k0x  =  1.57  cm-1  ,  k0y  =  2.39  cm-1  ,  ^kz0  -  0.732  cm"'  ,  and  Af  =  1.13x1 0-2  .  The 

solution  to  the  generalized  envelope  equations  shows  that  the  semiaxes  of  the  elliptical 
beam  remain  roughly  constant  with  small  oscillations,  that  the  orientation  of  the  ellipse 
twists  periodically  with  an  amplitude  of  ten  degrees,  and  that  the  normalized  rotation 
flow  velocities  ax  and  ay  oscillate  with  the  magnet  periodicity.  Shown  in  Figs.  3.5(a) 

and  3.5(b),  the  dotted  curves  are  the  envelopes  and  angle  of  the  beam  ellipse  obtained 
from  the  PFB2D  simulation.  It  is  worthwhile  pointing  out  that  the  normalized  velocities 
jux ,  /uy  ,  ax  and  ay  vanish  at  s  =  0  which  makes  it  a  natural  matching  point  for  a 
parallel-flow  beam  with  negligibly  small  emittance  (Bhatt  and  Chen,  2005). 
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Figure  3.5:  Plots  of  (a)  envelopes  <3(5)  and  b(s),  (b)  twist  angle  <?(s),  and  (c) 
normalized  rotational  velocities  b(s)ax(s)  and  a(s)ay  (s)  versus  the  axial  distance 

s  for  the  relativistic  twisted  ellipse-shaped  beam  in  Table  3.1,  Column  4.  The 
solid  and  dashed  curves  are  the  generalized  envelope  solution,  whereas  the  dotted 
curves  are  from  the  PFB2D  simulation. 
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Figure  3.6  Plots  of  10,000  particles  (a  sample  of  the  5  x  105  particles  in  the  PFB2D 
simulation)  in  the  ( x,y )  plane  and  ( x,dy/ds )  plane  for  five  snapshots  within  one 
period:  s/S  -  9.0, 9.25, 9.5, 9.75  and  10.0  for  the  same  beam  shown  in  Fig.  3.5. 


The  PFB2D  simulation  also  shows  that  the  transverse  beam  distribution  preserves  the 
equilibrium  profile  as  it  propagates.  In  Fig.  3.6,  10,000  particles  (a  sample  of  the 
5 xlO5 particles  in  the  PFB2D  simulation)  are  plotted  in  the  (x,^)  plane  and  (x,dy/ds) 
plane  for  five  snapshots  within  one  period:  s/S  =  9.0, 9.25, 9.5, 9.75  and  10.0  for  the  same 
beam  shown  in  Fig.  3.5.  This  also  suggests  that  the  beam  equilibrium  is  stable,  although 
we  have  not  done  a  full  stability  analysis  of  the  beam  equilibrium. 

These  high-intensity  twisted  elliptic  beam  equilibria  exist  over  a  wide  region  of 
parameters  ranging  from  the  nonrelativistic  to  relativistic  regimes.  In  addition  to  the 
relativistic  elliptic  beam  equilibrium  for  discussed  in  Figs.  2.3.5  and  2.3.6,  we  present 
three  additional  examples  in  Table  2.3.1  for  ribbon-beam  microwave  devices  which  have 
applications  ranging  from  wireless  communications,  to  high-power  radar  (e.g.,  AEGIS 
radar),  and  to  directed  energy. 

3.3.  Spin-off  of  Ribbon-Beam  Technology  for  Commercial,  Defense  and  Scientific 
Applications 

We  are  spinning  off  the  ribbon-beam  technology  to  Beam  Power  Technology,  Inc.,  a 
high-tech  start-up  company  whose  mission  is  to  apply  it  in  commercial,  defense,  and 
scientific  applications.  The  nonrelativistic  elliptic  beam  corresponds  to  a  beam  design  for 
a  high-efficiency  200  W  ribbon-beam  amplifier  (RBA)  which  is  being  developed  at 
Massachusetts  Institute  of  Technology  (MIT)  for  communication,  whereas  the  mildly 
relativistic  elliptic  beam  is  proposed  for  a  high-power,  high-efficiency  RBA  for  radar 
applications.  The  two  relativistic  elliptic  beams  are  proposed  for  International  Linear 
Collider  (ILC)  and  direct  energy  application,  respectively. 

4.  Chaotic  Particle  Motion  and  Halo  Formation  in  Small  Aperture  Systems  (Zhou 
and  Chen,  2005) 

It  is  important  to  identify  the  mechanisms  of  chaotic  particle  motion,  beam  halo 
formation,  and  beam  loss,  and  to  develop  techniques  to  control  them  in  particle  beam 
devices  in  general  and  HPM  devices  in  particular.  Under  the  auspices  of  the  present  grant, 
we  carried  out  a  parametric  study  of  chaotic  particle  motion,  beam  halo  formation,  and 
beam  loss  in  an  intense  charged-particle  beam  in  a  small-aperture  alternating-gradient 
focusing  system  (Zhou  and  Chen,  2005).  The  mechanism  for  chaotic  particle  motion, 
beam  halo  formation,  and/or  beam  loss  is  due  to  the  image  charge  effects  of  a  pulsating 
elliptic  beam  in  a  small-aperture  circular  conducting  pipe. 

Figure  4.1  shows  a  Poincare  surface-of  section  plot  for  the  trajectories  of  100  test 
particles  for  propagation  over  400  lattice  periods  in  the  phase  plane  (x,.y)  with 

normalized  pipe  radius  R  =  R/a  =  4.0,  where  a  is  the  maximum  value  of  the  semi 
major  axis  of  the  pulsating  elliptic  beam.  The  system  parameters  are  chosen  to  be  ^  =  0.5, 

KS/sx  =10.0 ,  Sx=£y=  S ,  and  crv  =  80°  .  In.  Fig.  4.1,  the  scattered  particles  are  chaotic. 

They  form  a  halo  which  extends  to  the  circular  conductor  wall,  permitting  particles  to 
strike  the  wall. 
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Figure  4.1:  Poincare  surface-of  section  plot  for  the  trajectories  of  100  test 
particles  for  propagation  over  400  lattice  periods  in  the  phase  plane  (x,y)  with 

the  system  parameters  R-  R/ a  =  4.0 ,  77  =  0.5,  KS/ sx  =  10.0,  £x—8y—S,  and 
<rv  =  80°  (from  Zhou  and  Chen,  2005). 


Figure  4.2:  Plot  of  the  range  of  the  transverse  energy  Wx  for  test  particles 
propagating  over  2000  lattice  periods  for  R  =  3.8,  77  =  0.5,  KS/ex=  10.0, 
sx  =  s  =  e ,  and  <rv  =  80°  .  The  115  test  particles  are  loaded  uniformly  at  5  =  0 
between  the  interval  0  <WX  <1.15  along  the  x  axis  in  phase  space  (from  Zhou 
and  Chen,  2005). 


-  *  * 


A  more  detailed  description  of  how  such  chaotic  particle  motion  arises  is  illustrated  in 
Fig.  4.2,  which  shows  a  plot  of  the  range  of  the  transverse  energy  Wx  for  test  particles 

propagating  over  2000  lattice  periods  for  R  =  3.8,  ^  =  0.5,  KS/sx  =10.0,  ex=-e  —  e, 
and  <rv  =  80°  .  The  1 1 5  test  particles  are  loaded  uniformly  at  s  -  0  between  the  interval 
0  <WX  <1.15  along  the  x'  axis  in  phase  space. 

In  designing  on  an  elliptic  electron  beam  for  HPM  device  applications,  it  is  important 
to  understand  image  charge  effects  and  determine  the  parameter  space  where  the  image 
charge  effects  are  negligibly  small.  This  will  be  an  important  area  of  future  research. 

5.  Research  Initiative  on  Relativistic  Electron  Diodes  and  Beams  for  HPM  Device 
Applications 

In  a  HPM  device,  a  relativistic  diode  is  used  to  generate  a  high-power  electron  beam 
which  interacts  with  rf  structure  to  generate  intense  electromagnetic  waves.  The 
performance  of  a  HPM  device,  in  general,  depends  strongly  on  how  well  the  relativistic 
electron  diode  and  subsequent  beam  system  are  design.  Despite  of  decades  of  theoretical 
and  experimental  research  on  HPM  devices,  our  understanding  of  relativistic  electron 
diodes  and  beams  is  limited.  To  date,  relativistic  gun  design  relies  on  computer 
simulations. 

A  major  problem  in  a  typical  circular  relativistic  electron  gun  is  that  the  beam  density 
fluctuates  as  it  propagates  in  the  focusing  magnetic  field.  The  density  fluctuations  cause 
undesirable  effects  such  as  beam  mismatch  (or  beam  scalloping),  emittance  growth,  and 
beam  halo  formation,  all  of  which  degrade  the  performance  of  HPM  sources  using  such 
beam. 

We  have  taken  an  initiative  to  develop  the  equilibrium  theory  of  axisymmetric 
( d/d6  =  0 )  relativistic  electron  diodes  and  beams  for  HPM  device  applications.  In 
particular,  we  have  formulated  the  relativistic  electron  beam  generation  and  propagation 
problem.  By  introducing  the  electron  flow  and  magnetic  flux  surfaces  x =  rB$  ar,d 
\j/  =  rAg,  respectively,  where  B  =  V  x  A  is  the  total  magnetic  field  and  A  is  the  vector 
potential.  The  energy  conservation  along  flow  lines  is  expressed  as 

ymc2  -e0  =  s(x),  (5.1) 

where  —e  is  the  electron  charge,  <j>  is  the  electrostatic  potential,  and  y  is  the  relativistic 
mass  factor.  The  azimuthal  angular  momentum  conservation  is  expressed  as 

rmyVe--y/  =  L{x).  (5.2) 

c 

It  can  be  shown  that  the  problem  of  relativistic  electron  diodes  and  beams  is  fully 
described  by  (Stowell  and  Chen,  2005) 


V2</>  =  \mn, 


(5-3) 
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(5.4) 


*  *•  - 


A*  A  V0 

A  \f/  =  4  nern  — , 
c 

T-^V^  fcxeJ-Vltoxe,))  -  -fvj-V*  _l!M— A  v*. Vr 

(4;rew)  m  4;re«r  cr 

(5.5) 

where  A*  =  r— — —  +  -^—~  and  V  = - —  Vxx^-  +  Vee0  .  Equations  (5.1)-(5.5) 

dr  r  dr  dz  4  nen  r 

provide  a  general  description  of  axisymmetric  relativistic  electron  diodes  and  beams. 

Solving  them  will  yield  various  equilibria  for  relativistic  electron  diodes  and  beams, 

which  will  be  our  future  investigations. 
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